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Abstract 

We have redone a recent two-loop computation of the critical mass 
for Wilson fermions in lattice QCD by evaluating Feynman integrals 
with the coordinate-space method. We present the results for differ- 
ent types of infrared regularization. We confirm both the previous 
numerical estimates and the power of the coordinate-space method 
whenever high accuracy is needed. 



1 Introduction 



The lattice formulation is at present the best tool for a nonperturbative study 
of QCD. Its most widely used implementation is due to Wilson |l|]. Here, the 
lattice fermion doubling problem is solved by the introduction of a formally 
irrelevant second-order term that gives a mass of the order of the cut-off to 
the doublers. This term explicitly breaks chiral invariance and therefore the 
limit of zero mass can be obtained only by a fine tuning, which is to say 
that the bare and the renormalized quark mass no longer vanish together. 
Nowadays, there are new formulations in which a form of chiral invariance 
is preserved at the cost of a nonlocal action (see, for example, the recent 
reviews [0, 0, Q and references therein). Even though this represents a very 
important progress from a conceptual point of view, it is not clear if the new 
formulations can become practical for numerical simulations. 

In the Wilson formulation, the breaking of chiral invariance gives rise to 
severe technical difficulties. Indeed, operators that have the correct chiral 
properties are obtained as sums of many different lattice bare operators mul- 
tiplied by suitable renormalization constants. In principle, such constants 
should be computed nonperturbatively. In practice, this is often difficult and 
thus one resorts to perturbation theory. 

Lattice perturbation theory is particularly complicated because of the 
loss of Lorentz invariance and at present a two-loop calculation represents a 
formidable task. For one-loop computations a completely general algebraic 
algorithm was introduced in Ref. ^ for the pure gauge theory and extended 
in Ref. |^ to include Wilson fermions. This algorithm allows to express 
every one-loop integral with gluon and Wilson-fermion propagators in terms 
of a small number of basic constants that can be computed with arbitrary 
precision 0. For two- loop calculations a powerful method was introduced 
by Liischer and Weisz |^ for the pure gauge theory, the so-called coordinate- 
space method. Its generalization to full QCD is again a nontrivial task. Some 
general ideas were presented in Ref. M and will be extensively illustrated 



elsewhere |]T0| 



Here, we present the first two-loop computation with the coordinate- 
space method in lattice QCD with fermions. The purpose is twofold. First, 
we wish to check that the additional complications that are present in the 
method in the presence of fermions do not spoil its performance. Second, 
we want to compare it with the alternative momentum-space method, in 
which one works directly in momentum space. For this purpose, we repeat 
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the calculation of the critical mass for Wilson fermions reported in Ref. |TT 



It turns out that the coordinate-space method is quite precise, providing 
numerical expressions for the Feynman integrals with a precision of 10~^° 
or higher. The momentum-space method also works well: although not as 
precise as the coordinate-space method, still the final results have a relative 
precision of 10~^. 

In this paper, we will consider Wilson fermions. The dressed inverse 
fermion propagator has the form 

S'~^(p,mij) =if +mB + Mw{p) - T.^{p,mB, go), (1) 

where, setting the lattice spacing equal to one, 

= sinp^ , (2) 

f = E(2sin^)' , (3) 

Mw{p) = "-ff. (4) 

The additive mass renormalization Stub is obtained by requiring S'~^(0, 5mB) - 
0, i.e. 

S^(0,5mB,5(o) = 5mB. (5) 
This equation can be solved in perturbation theory by expanding 

oo 

S^(0,mB,^?o) = E^?o"S('^^- (6) 

n=l 

In the following we shall compute S^^^ and S^^^ for = 1, gauge group 
SU{N), and Nf fermionic flavour species. We will work in the Feynman 
gauge. 



2 One-loop result 

In Ref. we already reported the analytic one-loop expression for the 
fermionic self-energy S^. The first computation for the Wilson action in 
Feynman gauge was given in Ref. [jl2| and it was subsequently corrected 



in Ref. |IT3|.0 Our result is expressed in terms of three purely bosonic 
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Zq 0.154933390231060214084837208 
Zi 0.107781313539874001343391550 
Fn 4.369225233874758 



Table 1: Numerical values of the three constants Zq, Zi^ and Fq. 



^(1,0) 


0.08539036359532067914 


^(1,- 


-1) 


0.46936331002699614475 


^(1,- 


-2) 


3.39456907367713000586 




-1) 


0.05188019503901136636 




-2) 


0.23874773756341478520 


m- 


-2) 


0.03447644143803223145 




-3) 


0.13202727122781293085 


^(3,- 


-4) 


0.75167199030295682254 






- 0.01849765846791657356 


Yi 




0.00376636333661866811 


Y2 




0.00265395729487879354 


Y^ 




0.00022751540615147107 



Table 2: Numerical values of the constants appearing in the fermionic inte- 
grals. 

constants — their numerical values are given in Table |1| — and of 12 numerical 
constants that appear in the presence of Wilson fermions, whose numerical 
values are given in Table ^. Except for Fq, we report here 20 digits, but 
many more can be computed, and indeed must be computed for the numeri- 
cal implementation of the coordinate-space method, see Ref. f^. In practice, 
we have computed all constants but Fq with 60-digit precision [|10|. These 
constants are defined as follows. We define the integrals 



l-n {27T)^DF{k,mf)PDB{k,mb)'^'' 

^Notice, however, that Eq. (3.15) m Ref. jl^ contains a misprint: the correct result is 
given in Eq. (10b) of Ref. (l|. 
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where 



DB{k,mb) 
DF{k,mf) = + Mw{kf + m). 

When one of the n's vanishes, we shall not write it. Then 



^0 



B{1)\ 



m(,=0 



while 



>3 



lim 

m— >0 



4^(i;i,i)U=o 

lim [167r2i3(2) +logm2 + 7^], 



^(2,0) + — ^(logm2 + 7^-Fc 



(8) 

(9) 

(10) 
(11) 

(12) 



•^(1,1; 1,1,1) 



lm=0 



- ^(l,l;l,l,l,l)L=o 
1 ^(1,2;1,1,1)L=„, 



(13) 



where we have set m = = mf,. Also, in Table ^the numerical values refer 
to the massless case m = 0. 
At one-loop order 



Ed) 



N 



i=l 



.(1) 



(14) 



where c\^^ are the contributions of the two diagrams illustrated in Fig. |1|. In 
terms of the basic integrals defined above, they are given by 



'-2 



f -^(1,0). 



The numerical values are reported in Table |^ and compared with the results 
of Ref . |TT[ , obtained by using the momentum-space method. The agreement 
is excellent. Summing up the two contributions we obtain 



^ cf ^ = -[Zo + 2 J^(l, 0)] ^ -0.16285705871085078618 . (15) 

1=1 



1 



2 



Figure 1: One- loop diagrams contributing to S^. Wavy (solid) lines repre- 
sent gluons (fermions). 



i 




1 


-0.15493339023106 




-0.15493339023106021408 


2 


-0.007923668480(2) 




-0.00792366847979057210 



Table 3: Coefficients c'p . For each of them we report in the ffist line the 



result of Ref. obtained by means of a momentum-space integration, 

and in the second line our result, obtained by means of the coordinate-space 
method. 



The constant is in excellent agreement with the result of Ref. [11 
Y^icf^ = 0.162857058711(2). 



3 Two-loop result 



At two loops there are 26 diagrams that have been drawn in Fig. 0. They 
are numbered as in Ref. ||ll| in order to simplify the comparison. The i-th 
diagram gives a contribution of the form 



A = (iv' - 1) 



• (2) J_ (2) 

^1,* + jY2 ^2,i 



+ 



N 



f M 



N 



(16) 



In Table ^ we report the results given in Ref. |TTj| and those obtained here 
by using the configuration-space method. When we have not reported an 
error, the precision we achieve is much higher than the reported digits. This 
occurs in general when the diagram is the product of one-loop integrals. All 
results are in agreement with what had been presented in Ref. [|ri|] . Only for 
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28 



Figure 2: Two-loop diagrams contributing to T,^. Wavy (solid, dotted) lines 
represent gluons (fermions, ghosts). Crosses denote vertices stemming from 
the measure part of the action; the solid circle is the one-loop fermion-mass 
counterterm. 
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I 


(2) 
"l.i 


2,i 


,(2) 
3,i 


3 


0.002000362950707492 


-0.0030005444260612375 







0.0020003629507074987148 


-0.0030005444260612480722 





4 


0.00040921361(1) 


-0.00061382041(2) 







0.0004092136068803147865 


-0.0006138204103204721798 





5 






















6 


-0.0000488891(8) 


0.000097778(2) 







-0.0000488853119(2) 


0.000097770623(5) 





7+8+9+10+11 


-0.013927(3) 
-0.01392647740 (2) 


0.014525(2) 

0.0145250053341618950704 






12+13 








0.00079263(8) 










0.000792647(2) 


14+15+16+17+18 


-0.005753(1) 
-0.00575248584(1) 


0.0058323(7) 
0.005832127004694453 






19+20 










0.000393556(7) 
0.000393556(4) 


21+22+23 


0.000096768(4) 
0.0000967648(2) 


-0.000096768(4) 
-0.0000967648(2) 






24 






















25 


0.00007762(1) 
0.000077613106(4) 


-0.00015524(3) 
-0.000155226212(8) 






26 


-0.00040000(5) 
-0.00039997586(1) 










27 





-0.000006522(1) 










-0.0000065203(1) 





28 


0.0000078482(5) 
0.0000078480652722033294 


-0.000015696(1) 
-0.0000156961305444066589 






Total 


-0.017537(3) 
-0.0175360218(2) 


0.016567(2) 
0.0165663304(2) 


0.00118618(8) 
0.001186203(6) 



Table 4: Coefficients c^i\, c^2i For each of them we report in the 



ffist hne the result of Ref. |]TT|, obtained by means of a momentum-space 
integration, and in the second line our result, obtained by means of the 
coordinate-space method. 
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diagram 6 there is apparently a (very) small underestimation of the error, 
which is negligible in the sum of all contributions. 

In Table ^ diagrams are grouped together in order to obtain infrared- 
convergent results. This is necessary for the numerical implementation of 
the momentum-space method. In our case, however, we have followed a 
different strategy considering an infrared regularization. This allows us to 
compute each Feynman diagram separately. We have used four different 
infrared regularizations: 

(a) We simply introduce a regulator in the propagators as in Eq. (||). Ex- 
plicitly, for the gluon (Ab(A;)) and for the fermion (Ap'(A;)) propagator 
we use: 

Ab(A;) = „ X y (17) 
Dp{k, m) 

(b) We regularize the gluon propagator as in Eq. (|13), but use instead the 
correct Wilson-fermion propagator 

^^^^^ - k^ + {M^{k)+my ■ ^^^^ 

(c) We regularize the Wilson fermion as in (a), but use the massless prop- 
agator for the gluon. 

(d) We regularize the Wilson fermion as in (b) and the massless propagator 
for the gluon. 

The result for each diagram will be indicated by D^^^ where (x) refers to the 
chosen infrared regularization. 

In Table ^ we report the result for each diagram computed in the regu- 
larization (a). As expected, the divergences cancel when the appropriate 
combinations of integrals are considered. The finite results are reported in 
Table |. 

Then, we consider regularization (b). In Table |^ we report, for two groups 
of graphs, the difference Df* —D^f' , where the superscripts refer to the chosen 
regularization. Note that, if we use the correct fermion propagator there are 
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i Divergent Part 



10 
11 



V 2567r 

( — 

V 38"- 



\ 256- 



3Zo ^ 

16^ ) 


1 log 


-0.01347769288130718 


0.0146260063341618960704 





^ log 




0.001000649213311037 








327r2 ; 


log 


-0.00347141693(2) 








- -£2^1 

327r^ ) 


1 log 


0.000020884775(4) 








\ log 




0.002001098426622074 









12 f^ + ii-ii^log^^ 

" V 327r2 167r2 / ^ 

13 i;;if^_iii-ii)iog,.= 



0.000735684385222637 
0.000056962(2) 



14 
15 
16 
17 
18 



\ 512 



( 256 



32 
7Zo 



) logr 

V 5127r2 647r2 y 
f ~S l logm^ 



-0.006494346068160373 

-0.00028296107(1) 

-0.000055474799728634 

0.0001912446811(2) 

-0.000110949699467269 



0.006832127004694453 











19 

20 



J[i.-i] 



-0.000040789541774414 
0.000434346(4) 



21 ( 




I f^logr„2 


0.001606284825541242 


-0.001606284825541242 





22 ( 


\ iV2 ; 




0.0006016205(2) 


-0.0006016206(2) 





23 I 




-0.002011040464066014 


0.002011040464066014 






Table 5: Results for the infrared regularization (a) at two loops. Graphs 
whose divergent contributions sum up to zero have been grouped together. 
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_Nj_ ( 1 _ _ J^(l.-l) J^(1.0) \ 

AT \^16m7r2 327r2 SmTr^ 4^2 j 

Nf_ ( 1 ^(1-1) _ _ J^(1,0) A 



22 


(1- 


1^) 


( -Za 


jr(l,0)\ 


23 


(1- 


^) 


( -Za 
\ 3271^ 


^(1:0) \ 
167r2 y 



Table 6: Difference D\ — D^. Graphs whose contributions sum up to zero 
are grouped together. 



21 
22 

23 




wil 



5Za 
1927r^ 

5Za 
1927r^ 



I 5J^(1.0) 
5J^(1 0) 



Table 7: Difference Df^ — Di' . Regularization (x) = (c) is obtained by 
taking Pwn = 0, regularization (x) = (d) by taking P^u = 1. 

additional divergences. As expected, the sum of the terms reported in Table | 
vanishes. 

As a last check, we have repeated the calculation of diagrams 21, 22, 
23 using the regularizations (c) and (d). In Table |^ we give the differences 
jji.^) — /^l") for X = c,d. The case x = c is obtained by setting P^n = 0, while 
X = d s obtained by setting P^u = 1. Again, we note that the individual 
diagrams are different, but the sum of the contributions vanishes. 

In conclusion, the coordinate-space method is a very efficient tool for 
the computation of (infrared- finite and infrared-divergent) two- loop Feynman 
integrals also in the presence of Wilson fermions. The method proposed in 
Ref. P really works. 
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